Velocity Renormalization and Carrier Lifetime in Graphene 
from Electron-Phonon Interaction 
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We present a first-principles investigation of tlie piionon-induced electron self-energy in graphene. 
The energy dependence of the self-energy reflects the peculiar linear bandstructure of graphene 
and deviates substantially from the usual metallic behavior. The effective band velocity of the 
Dirac fermions is found to be reduced by 4-8 %, depending on doping, by the interaction with 
lattice vibrations. Our results are consistent with the observed linear dependence of the electronic 
linewidth on the binding energy in photoemission spectra. 
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The recent fabrication of single-layer graphene [l[ has 
attracted considerable interest because low-energy charge 
carriers in this material have dispersion curves similar to 
Dirac fermions with zero rest mass and constant group 
velocity 0, [H . Because of the peculiar electronic struc- 
ture of graphene, electrons and holes exhibit exception- 
ally large mobilities, and the density of states can be 
tuned over a wide range by applying a gate voltage [1, [3] . 
These properties make graphene a promising candidate 
for new-generation electronic and spintronic devices. 

Angle-resolved photoemission spectroscopy (ARPES) 
is used as a powerful tool for investigating quasipar- 
ticle behavior with extremely fine energy and momen- 
tum resolution The photoelectron intensity pro- 
vides information about the energy vs. momentum dis- 
persions of the charge carriers and the associated life- 
times. Recent photoemission experiments performed on 
graphene showed a peculiar dependence of the hole life- 
time on the binding energy, as well as a significant ve- 
locity renormalization [5] . The measured carrier lifetime 
has been discussed within a model including three differ- 
ent decay channels: electron-phonon (e-ph) scattering, 
electron-plasmon scattering, and electron-hole pair gen- 
eration [5]. The linear dependence of the linewidth on 
the binding energy was attributed to the generation of 
electron-hole pairs. The phonon-induced lifetime was as- 
sumed to be energy-independent as found in conventional 
metallic systems [5|. A subsequent theoretical work an- 
alyzed the carrier lifetimes in graphene by adopting a 
two-dimensional electron-gas model, and concluded that 
the experimental results could be explained without in- 
voking the e-ph interaction . 

In this work we investigate the e-ph interaction in 
graphene within a first-principles approach. We calcu- 
lated the electron self-energy arising from the e-ph inter- 
action using a dense sampling of the scattering processes 
in momentum space, and we extracted the velocity renor- 
malization and the carrier lifetimes from the correspond- 
ing real and imaginary parts of the self-energy, respec- 
tively. Our analysis shows that the self-energy associ- 
ated with the e-ph interaction in graphene is qualitatively 
different from that found in conventional metals. The 



imaginary part of the self-energy shows a linear energy 
dependence above the phonon emission threshold, which 
directly reflects the bandstructure of graphene. The real 
part of the self-energy leads to a Fermi velocity renormal- 
ization of 4-8 % depending on doping. We further pro- 
pose a simple analytical model of the electron self-energy 
capturing the main features of our first-principles calcula- 
tions. Our calculation allows us to assign the low-energy 
kink in the measured photoemission spectrum and part of 
the linear energy-dependence of the electronic linewidths 
to the e-ph interaction. 

The e-ph interaction in graphene is treated within the 
Migdal approximation [7]. The contribution to the elec- 
tron self-energy S„k(-E'; T) arising from the e-ph interac- 
tion at the temperature T is 0, S] : 
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where e„k is the energy of an electronic state with band 
index n and wavevector k, and hLOc^u the energy of a 
phonon with wavevector q and branch index v. fnk and 
nqjy, are the Fermi-Dirac and Bose-Einstein distribution 
functions, respectively. The integration extends over the 
Brillouin zone (BZ) of graphene of area Abz and the 
sum runs over both occupied and empty electronic states 
and all phonon branches. The e-ph matrix element is 
defined by g,„„,,,(k, q) = (mk -|- q|AVqi.|nk), AVq^ be- 
ing the change in the self-consistent potential due to a 
phonon with wavevector q and branch index v, while 
|nk), |mk-f q) indicate Bloch eigenstates. Equation ([1]) 
takes into account the anisotropy of the e-ph interac- 
tion in k-space, as well as retardation effects through the 
phonon frequency in the denominators. 

The electronic structure was described within the local 
density approximation to density-functional theory 
Valence electronic wavefunctions were expanded in a 
plane-waves basis [13] with a kinetic energy cutoff of 
60 Ry. The core-valence interaction was treated by 
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means of norm-conserving pseudopotentials flU . Lattice- 
dynamical properties were computed through density- 
functional perturbation theory ■ We modeled an iso- 
lated graphene by a honeycomb lattice of carbon atoms 
within a periodic supercell. The graphene layers were 
separated by 8.0 A of vacuum [l^, and the relaxed C- 
C bond-length was 1.405 A. Doped graphene was mod- 
eled by varying the electronic density and introducing 
a neutralizing background charge. We first calculated 
electronic and vibrational states and the associated e- 
ph matrix elements on 72 x 72 k-points and 12 x 12 q- 
points in the BZ of graphene. Then, we determined 
the quantities needed to evaluate the self-energy given 
by Eq. ([1]) on a significantly finer grid of 1000 x 1000 k 
and k-l-q points in the irreducible wedge of the BZ by us- 
ing a first-principles interpolation based on electron and 
phonon Wannier functions 14, [l^ [3]. The fine sam- 



pling of the BZ was found to be crucial for convergence 
of the self-energy. In the calculation of the self-energy we 
used a broadening parameter S of 10 meV, comparable 
with the resolution of state-of-the-art photoemission ex- 
periments fSl . The calculations were performed with the 
electron and phonon occupations [Eq. ([T])] corresponding 
to T = 20 K to make connection with the ARPES ex- 
periment [5j. In what follows, we discuss the computed 
electron self-energy by focusing on a straight segment 
perpendicular to the FK direction and centered at the K 
point in the BZ [Fig.[T]. 

We note here that within 2.5 eV from the Dirac point, 
the angular dependence of the self-energy is insignificant 
(at fixed energy E) [l3|. As a consequence, the e-ph 
couphng parameter A„(k) — ~ d KeT^n'k.iE) / dE\^^^ is 
isotropic in k-space. 

Figure [T] shows the calculated imaginary part of the 
self-energy (which is closely related to the linewidth) as 
a function of carrier energy, corresponding to three repre- 
sentative situations: intrinsic, electron-doped, and hole- 
doped graphene. We here considered doping levels cor- 
responding to 4-10^^ cm^^ electrons or holes. The corre- 
sponding Fermi levels were found to lie at -1-0.64 eV and 
—0.66 eV from the Dirac point, respectively. In the in- 
trinsic system, we found that the electron linewidth due 
to e-ph interaction is negligible (< 1 meV) within an en- 
ergy threshold ?iWph for the emission of optical phonons 
{fiLOp\y K, 0.2 eV being a characteristic optical phonon fre- 
quency) , while it increases linearly beyond this threshold 
[Fig.[lja)]. The scattering rate for electrons with energy 
below the optical phonon emission threshold is negligi- 
ble because (i) only optical phonons are effective in e-ph 
scattering and (ii) Pauli's exclusion principle blocks tran- 
sitions into occupied states. On the other hand, the lin- 
ear increase of the linewidth above the optical phonon 
energy relates to the phase-space availability for elec- 
tronic transitions, and reflects the linear variation of the 
density of states around the Dirac point in graphene. 
The energy dependence of the electron linewidths in the 
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FIG. 1: Calculated imaginary part of the electron self-energy 
arising from the e-ph interaction at T = 20 K (solid lines), for 
(a) intrinsic, (b) electron-doped, and (c) hole-doped graphene. 
The self-energy Ek(ek) was evaluated along the reciprocal 
space line segment shown in the upper-left corner. The Fermi 
level and the Dirac point are shown schematically in each 
case. We also show for comparison the imaginary part of the 
self-energy for a conventional metal (dashed lines) (Ref. IT^. 



electron-doped and the hole-doped systems [Fig.[TJb) and 
Fig. [He), respectively] can be rationalized by a similar 
phase-space argument. We denote by E'd the energy of 
the Dirac point with respect to the Fermi level. For def- 
initeness, we here consider the electron-doped situation 
(i?D < 0). When the energy of the hole is exactly equal 
to IE^dI + ?JWph (i.e., at — |i?D| — in Fig. [T|), there are 
no allowed final states for electronic transitions through 
optical phonon emission, resulting in a vanishing scat- 
tering rate at zero temperature. As the hole energy de- 
parts from I -Ed I + ?iWph, the linewidth increases linearly, 
and exhibits a characteristic dip at the Fermi level. The 
latter feature corresponds to forbidden phonon emission 
processes, as discussed above for intrinsic graphene. The 
calculated energy dependence of the electron linewidth 
deviates substantially from the standard result which ap- 
plies to conventional metals (Fig. [U dashed line) [l8|. 
The latter consists of a constant scattering rate above 
the phonon emission threshold, and fails in reproducing 
the features revealed by our ah initio calculations. 

Figured] shows the real part of the electron self-energy 
arising from the e-ph interaction, for intrinsic and for 
electron-doped graphene. The behavior of the hole- 
doped system is qualitatively similar to the electron- 
doped case. While in conventional metals the real part of 
the self-energy decays at large hole energies {E < — TiWph) 
[Fig.[2fb), dashed line], the self-energy in graphene shows 
a monotonic increase in the same energy range [Fig.[2l^b), 
solid line] . Since the wavevector dependence of the self- 
energy in graphene within a few eV from the Fermi level 
is neghgible [i.e., S„k(-E) ~ S„(E)] [13], we obtained the 
quasiparticle strength Z„k = (1 — 9 ReI]„k/9E)~^ by 
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FIG. 2; Calculated real part of the electron self-energy arising 
from the e-ph interaction at T = 20 K (solid lines), for (a) 
intrinsic and (b) electron-doped graphene. The self-energy 
was evaluated along the reciprocal space line segment shown 
in Fig. [1] The corresponding velocity renormalization (w„k — 
u°i^)/w„k is shown in panels (c) and (d), respectively. We also 
report, for comparison, the real part of the self-energy and 
the velocity renormalization for a conventional metal (dashed 
lines) (Ref. [H). At variance with conventional metals, the 
group velocity in graphene shows additional dips when the 
carrier energy is |i?D| +i^ph (arrows), reflecting the vanishing 
density of states at the Dirac point. 



evaluating {1 — d ReS]„k(ek)/c^ek;)~^- In all cases con- 
sidered, the e-ph interaction was found to reduce the 
non-interacting quasiparticle strength down to at most 
ZnU = 0.93 at the Fermi level. This suggests that a 
quasiparticle picture is still appropriate at low energy, the 
e-ph interaction largely preserving the weakly perturbed 
Fermi-liquid behavior. The quasiparticle strength is re- 
lated to the velocity renormalization through 1 — = 
{vnk - v^■^J/vnk, and w„k being the non- interacting 
and the interacting velocity, respectively. The velocity 
renormalization is plotted in Fig. ^c) and Fig. [IJd) for 
the intrinsic and the electron-doped system, respectively. 
The velocity renormalization at the Fermi level was found 
to increase with the doping level, and amounts to —4 %, 
—8 %, and —6 % in the intrinsic, the electron-doped, and 
in the hole-doped system considered here. Our results in- 
dicate that the velocity of Dirac fermions in graphene is 
affected by the e-ph interaction. This bears important 
implications for the transport properties of graphene- 
based electronic devices. 

In order to provide a simplified picture of the e-ph in- 
teraction in graphene, we analyzed the various e-ph scat- 
tering processes contributing to the electron lifetimes. 
We repeated our calculations by restricting either the 
energy huqi, or the momentum transfer q in Eq. ([T]) to 
limited ranges. When only the in-plane optical phonon 
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FIG. 3: Comparison between the electron self-energy obtained 
from a first-principles calculation (solid lines) and a single- 
parameter model (dashed lines): imaginary (upper panels) 
and real (lower panels) part for the intrinsic system (left) as 
well as for the electron-doped system (right). Note that the 
horizontal energy ranges differ from those shown in Fig. [T] 



modes between 174 and 204 meV are taken into account 
in Eq. ([T]) , the electron linewidth is found to deviate from 
the full ah initio result by 15% at most. In contrast, 
when the momentum integration in Eq. ([1]) is restricted 
to small regions around the high-symmetry points F and 
K, the linewidth is found to deviate significantly from the 
full calculation, indicating that a proper account of the 
entire BZ is essential. Based on this analysis, we devised 
a simplified single-parameter model of the e-ph interac- 
tion in graphene. We assumed: (i) linear electronic dis- 
persions up to a few eV away from the Dirac point [l^ . 
(ii) an Einstein model with the effective phonon energy 
Titjph set to that of the highest degenerate zone-center 
mode, (iii) an effective e-ph vertex q, independent of the 
electron and phonon momenta [l7|. The e-ph matrix 
element g represents a free parameter in our simplified 
model, and has been determined by matching the model 
self-energy with the full ah initio result. Within these 
assumptions, and with the Fermi level set to zero, the 
imaginary part of the self-energy reads 



Im Y.{E) = 



~T6 



-a: 



2 „2 



\E- sgn{E) htOph - 



(2) 



whenever \E\ exceeds the characteristic phonon energy 
tiujpii, and vanishes otherwise. In Eq. ([2|), a is the lat- 
tice parameter in Bohr units, Uq = /fiv^ = 2.53 is 
the effective fine structure constant of graphene, and g 
is the average e-ph matrix element in Rydberg units. 
The fitting to our calculated ah initio self-energy gave 
g = 3.5 • 10^^ Ry. The real part of the model self-energy 
can be straightforwardly obtained from Eq. ([2]) through 
Kramers-Kronig relations. Figure [3] shows that this sim- 
plified model is in fairly good agreement with the full 
first-principles calculation. Therefore, despite its simplic- 
ity, our single-parameter model captures the qualitative 
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NERSC. Part of the calculations were performed using 
the Quantum-Espresso f20j and Wannier [2l| packages. 

Note added: After submission of this manuscript we be- 
came aware of a related work whose results are in agree- 
ment with our conclusions 22l |. 
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FIG. 4: Calculated width of the ARPES momentum distribu- 
tion curve for electron doped graphene (solid line) compared 
to the experimental result of Ref. 5: (dashed line). In our 
calculation, the Fermi level was set in order to simulate the 
sample with 2.1 ■ 10^"^ electrons/cm^ in Fig. 3 of Ref. 



features of the e-ph interaction in graphene. 

In Fig. 3] we compare our first-principles calculations 
with the width of the momentum distribution curve 
(MDC) measured by ARPES experiments at 20 K on 
graphene with a similar doping 5]. The width AA:„k of 
the MDC was calculated taking into account renormal- 
ization effects through Afc„k = .^nk 2Im S„k/^w„k 7|. 
Figure m shows that, contrary to previous findings 0|, 
the e-ph interaction plays a significant role in reducing 
the carrier lifetime in graphene, as it accounts for about 
a third of the measured linewidth at large binding en- 
ergies. The e-ph contribution to the width of the MDC 
is found to increase linearly at large binding energy, in 
agreement with experiment. 

In conclusion, we have computed from first-principles 
the velocity renormalization and the carrier lifetimes 
in graphene arising from the e-ph interaction and we 
have reproduced these results with a simplified model. 
The calculated energy-dependence of the phonon-induced 
electronic linewidths is shown to originate from the linear 
electronic dispersions. The renormalization of the Fermi 
velocity was found to be —4 % for intrinsic graphene and 
—8 % for an electron doping of 4 • 10^"^ cm^^, and is ex- 
pected to affect the mobility of graphene-based electronic 
and spintronic devices. 
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